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Abstract 

Background: In fabrication of ZnO-based low voltage varistor, Bi203 and Ti02 have been used as fornner and grain 
growth enhancer factors respectively. Therefore, the nnolar ratio of the factors is quit important in the fabrication. In 
this paper, modeling and optimization of Bi203 and Ti02 was carried out by response surface methodology to 
achieve maximized electrical properties. The fabrication was planned by central composite design using two 
variables and one response. To obtain actual responses, the design was performed in laboratory by the 
conventional methods of ceramics fabrication. The actual responses were fitted into a valid second order algebraic 
polynomial equation. Then the quadratic model was suggested by response surface methodology. The model was 
validated by analysis of variance which provided several evidences such as high F-value (153.6), very low P-value 
(<0.0001), adjusted R-squared (0.985) and predicted R-squared (0.947). Moreover, the lack of fit was not significant 
which means the model was significant. 

Results: The model tracked the optimum of the additives in the design by using three dimension surface plots. In 
the optimum condition, the molars ratio of Bi203 and Ti02 were obtained in a surface area around 1.25 point that 
maximized the nonlinear coefficient around 20 point. Moreover, the model predicted the optimum amount of the 
additives in desirable condition. In this case, the condition included minimum standard error (0.35) and maximum 
nonlinearity (20.03), while molar ratio of Bi203 (1.24 mol%) and Ti02 (1.27 mol%) was in range. The condition as a 
solution was tested by further experiments for confirmation. As the experimental results showed, the obtained 
value of the non-linearity, 21.6, was quite close to the predicted model. 

Conclusion: Response surface methodology has been successful for modeling and optimizing the additives such as 
Bi203 and Ti02 of ZnO-based low voltage varistor to achieve maximized non-linearity properties. 

Keywords: Optimization, ZnO-varistor, Modeling, RSM, Bi203, Ti02 



Background However, ZnO-based varistor has been demanded along 

Varistors are nonlinear electro-devices with a ceramics with the development of very-large- scale integration 

microstructure that are used as protectors in distribution electronics because it exhibits high nonlinear current- 

and energy transmission lines against voltage surge [1]. voltage (I-V) characteristics in lower voltage ranges 

In the past four decades, varistors based on ZnO and [4,5]. The non-linearity is expressed by I = K V" where K 

Sn02 have attracted attention because of their excellent is a constant, and a' is nonlinear coefficient (alpha) [6]. 

non-ohmic behavior and low leakage currents [2,3]. The alpha originates from microstructure of the varistor 

ceramics which is composed of conductive n-type ZnO 

grains and small amount of a few metal oxide additives 
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contain of Bi-rich intergranular, metal oxide and se- 
condary spinel phase and strictly influences on the 
alpha [4,6-11]. The role of Bi203, as a former, is quite 
important since it provides the medium for liquid- 
phase sintering, enhances the growth of ZnO grains, 
and finally stables the nonlinear current-voltage cha- 
racteristics of the varistor [12]. High sintering tem- 
perature is necessary for ZnO grain growth despite the 
fact that at this condition Bi203 tends to evaporate 
[13]. The melting point of Bi203 is 825°C, and the eu- 
tectic temperature of ZnO-Bi203 is only 740°C, thus a 
liquid phase is formed in the ZnO-Bi203 specimens 
below 800°C. As soon as the eutectic liquid is formed, 
the mass loss starts to increase which indicates the 
vaporization of Bi203 [14]. The sharp lost weight was 
reported above 1100°C since there was no reported 
peaks of |3-Bi203 at 1300°C [15,16]. On the other hand, 
Ti02 increases reactivity of the Bi203-rich liquid phase 
with the solid ZnO during sintering process which pre- 
vents Bi203 vaporization [13,17-19]. The phase equili- 
brium and the temperature of Uquid-phase formation 
are defined by the Ti02/Bi203 ratio [20]. According to 
the reports, the effect of Ti02 depends on Bi203 that 
means the additives interact in low-voltage varistor ce- 
ramics fabrication. To determine the effect of the inter- 
actions on the varistors' electrical properties, the molar 
ratios of the additives must simultaneously be consi- 
dered. To the best of our knowledge, there is no study 
on the interactions which optimize the ratio of the 
additives and maximize the alpha. Recently, response 
surface methodology (RSM) has been accepted for 
modeUng and optimizing of input intractable variables 
to achieve maximum yield product as output for pro- 
ductive process [21]. RSM is known as a semi-empirical 



method because the process could be optimized by 
using experimental results, and a group of mathema- 
tical and statistical techniques [22]. In this work, RSM 
was used for modeling and optimizing of molar ratio of 
Bi203 and Ti02 as additives to achieve the maximum 
value of the alpha for low voltage varistor. The experi- 
ments were designed by central composite design (CCD) 
to obtain the empirical results (actual). The results were 
used for regression and fitting process to fine an appro- 
priate model. The model was verified by several statistical 
techniques such as residual analysis, scaling residuals and 
prediction error sum of squares (PRESS). The model 
optimized the input additives and then maximized the 
alpha as output. In addition, the model predicted the 
desirable condition including minimum standard error 
and the maximum alpha which are validated by further 
experiments. The predicted samples were characterized by 
X-ray diffractometer (XRD), scanning electron microscope 
(SEM), variable pressure scanning electron microscope 
(VPSEM) and Energy-dispersive X-ray (EDX). 

Experimental 

Materials and methods 

The commercial chemical, ZnO (99.99%), Bi203 
(99.975%), Ti02 (99.9%), Sb203 (99.6%), Mn304 (98%), 
C03O4 (99.7%) and A1(N03)3 (100% ±2), were pro- 
vided from Alfa Aesar as starting powders. The pow- 
ders were weighed according to the experimental 
design of molar ratios (Table 1). The molar ratios of the 
powders were mixed, grounded in dry form and then ball 
milled in acetone for 24 h. During ball milling, agglo- 
meration was controlled by Zirconium oxide balls. After 
drying in hot oven for 8 h, the mixed powders were 
grounded and pressed into pellet forms of 10 mm in 



Table 1 Experimental-design contain of the actual variables, and actual response and model predicted values of the 
alpha 



Run 


ZnO 


Ti02 


BizOs 


CO3O4 


IVIn203 


Sb203 


AI(N03)3 


Alpha (Actual) 


Alpha (Predicted) 


1 


96.50 


1 


1 


0.5 


0.5 


0.5 


0.00094 


9.3 


10.1 


2 


96.00 


1.5 


1 


0.5 


0.5 


0.5 


0.00094 


3.9 


3.9 


3 


96.00 


1 


1.5 


0.5 


0.5 


0.5 


0.00094 


6.4 


74 


4 


95.50 


1.5 


1.5 


0.5 


0.5 


0.5 


0.00094 


10.5 


10.6 


5 


96.35 


0.896 


1.25 


0.5 


0.5 


0.5 


0.00094 


9 


7.9 


6 


95.65 


1.604 


1.25 


0.5 


0.5 


0.5 


0.00094 


5.6 


5.8 


7 


96.35 


1.25 


0.896 


0.5 


0.5 


0.5 


0.00094 


8.2 


7.7 


8 


95.65 


1.25 


1.604 


0.5 


0.5 


0.5 


0.00094 


11 


10.5 


9 


96.00 


1.25 


1.25 


0.5 


0.5 


0.5 


0.00094 


20 


20 


10 


96.00 


1.25 


1.25 


0.5 


0.5 


0.5 


0.00094 


19.6 


20 


11 


96.00 


1.25 


1.25 


0.5 


0.5 


0.5 


0.00094 


20.2 


20 


12 


96.00 


1.25 


1.25 


0.5 


0.5 


0.5 


0.00094 


20.7 


20 


13 


96.00 


1.25 


1.25 


0.5 


0.5 


0.5 


0.00094 


194 


20 
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diameter and 070 mm thickness at 200 MPa by a uniaxial 
presser machine. The disks were sintered in a box furnace 
(CMTS Model HTS 1400) for holding time of 2 h at 
1260°C. The heating and cooling rate were 5°C/min [23]. 
To determine DC current-voltage (I-V), both of the 
sintered sample surfaces were coated by silver electrodes. 
The I-V of the samples was measured by Keithley 236 
source meter. The samples were scanned with dc voltage 
from 0 to 100 V in step size of 2.5 V. The alpha was calcu- 
lated at Ji = 0.1 and J2 = 1 mA/cm^ by equation (1) as 
actual responses [6]. 

a = (logl2-logIi)/(logV2-logVi) (1) 

The breakdown voltage (Eb) was determined by mea- 
suring E at J = 0.75 mA/cm^ and the leakage current (Jl) 
was determined evaluating J at 0.8Eb where J (mA/cm^) 
is the current density and E is the electrical field (V/mm). 
To characterize the microstructure, the both surfaces of 
samples were polished by aluminum oxide powder. Then, 
they were etched at 160°C under sintering time with 
heating and cooling rate, 10°C/min. Phase analysis was 
conducted using XRD (PANalytical, Philips-Xpert Pro 
PW3040/60) with CuKa source. The sample were radiated 
with Ni-filtered CuKa radiation (\ = 1.5428) within the 26 
scan range of 20-80°. Surface morphology and elemental 
analyses of sintered samples were studied under SEM 
(JEOL JSM 6400) and VPSEM (LEO 1455) which attached 
to EDX. The samples was mounted on Al stub using car- 
bon paint and coated by gold layer. Average grains size of 
the ZnO in the varistor was evaluated by measuring 100 
grains in SEMs images. 

Experimental design 

The experimental design was carried out by CCD that 
used Design-Expert software version 8.0.7.1, Stat-Ease 
Inc., USA [24-26]. CCD is well fixed for fitting a quad- 
ratic surface that usually works well for optimization 
process [27-29]. The variables number, level of variables 
and number of responses are determined as input of ex- 
perimental design. In this case, the molar ratio of Bi203 
and Ti02 was selected as effective variables in vicinity of 
their optimum while alpha was the response as output. 
The CCD transformed the variables and the response to 
the terms of code values (Table 2) because the units and 



Table 2 The variables and employed levels in the CCD for 
ZnO low voltage varistor fabrication 









Level of variables 




Symbol 


Unit (%) 


Low (-1) 


Middle (0) 


High (+1) 


Bi203 (xi) 


mol 


1.0 


1.25 


1.5 


Ti02 (X2) 


mol 


1.0 


1.25 


1.5 



range of variables were different. The coded values 
spaces are ±1 from the center (0.0) and the star points 
are usually located 'a distance from the center [29,30]. 
In the design, there N experiment that includes the fac- 
torial points (2^), the axial points (2n), and the center 
points Co or replications as the equation N = 2^ + 2n + 
Co which is 13, 4, 4 and 5 respectively. The replications 
are used to measure experimental error [31]. As a result, 
the experimental design is presented as a design matrix 
with 'n' column and 'N' row in Table 1. Where, each co- 
lumn corresponds to a particular variable, e.g. Xi and X2 
which arranged in order to increase factor number from 
left to right. The rows are experiments runs because 
each one contains the descriptions of an experiment. 
Additionally, the design constructs a matrix of actual 
responses that obtained by the experiments. 

The RSM description 

The RSM develops an adequate functional relationship 
between input variables and interested responses by low- 
degree approximation of the polynomial models such as 
the second degree model (Eq. 2) [32]. 

i=l i=l i=l ;■=/+! 

(2) 

where Y is the interested response, is a constant term, 
/3i is the coefficient of the linear terms, represents the 
coefficient of the quadratic terms, is the coefficient of 
the interaction terms while Xi are control variables and 
's is a random experimental error [31]. For system with 
two factors, the model is described by equation (3), 

(3) 

where is the experimental single response, Xi and X2 
are the coded factors (Table 2), is the intercept term, 
Pi and ^2 are slopes with respect to each of the two fac- 
tors, and P22 are curvature terms, and ^12 is the 
interaction term. To estimate the y^s, the fitting 
process provides the sufficient data by regression 
tools [33,34]. In the process, the actual responses are 
fitted to the polynomial models by sequential model 
sum of squares (SMSS) [33,34]. The SMSS compares 
the linear, two-factor interaction (2FI), quadratic and 
cubic models by using the statistical significance of 
adding new model terms, step by step in increasing 
order [35]. To select the provisional model, the lack 
of fit of those models is compared by minimum p- 
values and PRESS. The other assessments to select 



Abdollahi et al. Chemistry Central Journal 2013, 7:137 
http://journal.chennistrycentral.conn/content/7/1/137 



Page 4 of 10 



the best provisional model are maximum adjusted R- 
squared (RAdj) and predicted R-squared (Rpred) [36]. 
The p-value is one of the most important evidences 
which was used to study significant effect of the para- 
meters [33]. In addition, the lack-of-fit test diagnoses 
how well each terms of the full model fit the data that 
pillared by statistical parameters such as RAdj> Rpred and 
PRESS [36,37]. Therefore, the provisional model with 
minimum p-value and PRESS and also maximum RAdj> 
Rpred is selected to investigate in details. The details are 
provided by using analysis of variance (ANOVA) which 
contains a collection of terms statistical evidences. The 
ANOVA determines the significance of intercept, linear, 
interaction and square terms of the provisional model 
by using minimum p-value. For more evaluation, the 
normality of residuals, constant error and residual out- 
lier is checked by various diagnostic plots [38]. The vali- 
dated model, the relationship between variables and 
response, is created in coded and actual variables. The 
model indicates the effect of linear, quadratic and the 
parameters interactions on the interested response. The 
effects are presented by estimated coefficients and the 
related positive and negative signs (+, -). The coeffi- 
cients are specific weight of the parameters in the model 
while the signs (+) and (-) operate as synergistic and 
antagonistic effects on the response [39]. Then the 
optimization process investigates combination of va- 
riables levels that produces the maximum response to a 
surface area simultaneously. Moreover, the model pre- 
dicts the yield of product in specific condition such as 
individual standard error, the range of variables and 
responses. The prediction could be performed by fur- 
ther experiments. 



Results and discussion 

Modeling 

In the fitting process, the residuals are produced from 
difference between actual and predicted values. Standard 
error is a great tool to determine the residuals outlier 
and also the scope of models prediction [40]. Figure 1 
indicates the standard error contour plot of the experi- 
ment design which displays Bi203 versus Ti02 molar ra- 
tio. The bright area has relatively low standard error that 
interested for the modeling and optimization process. 
However, the darker shading corners represent higher 
standard errors which are dangerous to extrapolation. 

Table 3 indicates that the SMSS compared four models 
to recommend a proper model [35,36]. As shown, the 
quadratic model was suggested as the best provisional 
model. The suggestion based on the lowest standard 
deviation, p-value and PRESS and the highest RAdj? Rpred 
values. Moreover, the RAdj (0.991) was in reasonable 
agreement with (<0.20) Rpred (0.985) which confirmed 
the model sufficiency. Therefore, the authors selected 
the suggested model to investigate in details by using 
ANOVA. 

Table 4 shows the ANOVA of the provisional model 
which included the useful statistical evidences about the 
terms {xi^ and xl) in details. As shown, the 

prob > F of the terms was less than 0.05 which con- 
firmed the high significance of the terms. Moreover, 
the models F-value was 153.6 that indicated great 
significance for the model. In addition, the very low 
value of the model p-value confirmed the significance. 
Futhermore, R-squared (R) provides a measure of 
how much variability in the observed response values 
can be explained by the experimental factors and their 
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Table 3 The sequential model fitting summary for the actual responses which shows statistics conformation of the 
regression process, DF is degree of freedom 



Source 




Sum of squares 


DF 


Mean square 


F-Value 


p-value 


Remark 


Sequential model sum of squares 


Mean vs Total 




2064.18 


1 


2064.18 








Linear vs Mean 




12.07 


2 


6.04 


0.13 


0.8818 




2FI vs Linear 




22.31 


1 


22.31 


0.44 


0.5216 




Quadratic vs 2FI 




447.1 1 


2 


223.55 


356.58 


< 0.0001 


Suggested 


Cubic vs Quadratic 




1.56 


2 


0.78 


1.38 


0.3325 


Aliased 


Residual 




2.83 


5 


0.57 








Total 




2550.06 


13 


196.16 








Source 




Sum of squares 


DF 


Mean square 


F-Value 


p-value 


Remark 


Lack of fit tests 


Linear 




472.78 


6 


78.80 


307.73 


< 0.0001 




2FI 




450.47 


5 


90.09 


351.85 


< 0.0001 




Quadratic 




3.36 


3 


1.12 


4.38 


0.0938 


Suggested 


Cubic 




1.80 


1 


1.80 


7.03 


0.0569 


Aliased 


Pure Error 




1.02 


4 


0.26 








Source 




Std.Dev. 




Rpred 


R 


PRESS 


Remark 


Model summary statistics 


Linear 




6.88 


0.025 


-0.170 


-0.571 


763.55 




2FI 




7.08 


0.071 


-0.239 


-1.195 


1066.45 




Quadratic 




0.79 


0.991 


0.985 


0.947 


25.53 


Suggested 


Cubic 




0.75 


0.994 


0.986 


0.760 


116.85 


Aliased 


interactions. In this study, the R (0.991) indicated that the 
model was capable of accounting for more than 99.1% of 
the variability in the responses. The RAdj (0.985) was in 
reasonable agreement with (<0.20) the Rpred (0.947) which 
confirmed the aptness of the model. The pure errors such 
as experimental errors were minimal as the lack of fit 
(0.094) was not significant or the model was fit well. 


Therefore, ANOVA confirmed the adequacy of the quad- 
ratic model that could be used to navigate the design 
space. 

The normality of residuals, constant error and residual 
outliers were checked by various diagnostic plots. The 
normal probability plots of the studentized residuals as 
one of the most important diagnostic plots, was provided 


Table 4 Analysis of variance (ANOVA) for response surface quadratic model, MS is mean Square, DF is degree of 
freedom and SS is sum of squares while and X2 introduce in Table 2 


Source 


SS 


DF 


MS 


F-Value 


p-value (Prob > F) 




Suggestion 


Model 


481.5 


5 


96.3 


153.6 


< 0.0001 




significant 


Xl 


4.6 


1 


4.6 


7.4 


0.0299 






^2 


7.4 


1 


7.4 


11.9 


0.0108 






X1X2 


22.3 


1 


22.3 


35.6 


0.0006 






X? 


298.7 


1 


298.7 


476.4 


< 0.0001 






xi 


205.5 


1 


205.5 


327.7 


< 0.0001 






Residual 


4.4 


7 


0.6 










Lack of Fit 


3.4 


3 


1.1 


4.4 


0.0938 




not significant 


Pure Error 


1.0 


4 


0.3 




< 0.0001 






Cor Total 


485.9 


12 






0.0299 






Std.Dev. 




Rpred 




PRESS 


c.v.% 




Adeq Precision 


0.792 


0.985 


0.947 


0.991 


25.5 


6.284 




29.9 
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by software default (Figure 2). The plot presents per- 
centage of normal% probability versus internal studentized 
residuals. The studentized residual is an important tech- 
nique in the detection of residuals outliers in regressions 
[41]. As the plot demonstrates, the residuals followed a 
normal distribution that implies the points follow a 
straight line. 

As another diagnosis, Figure 3 illustrates the predicted 
response values versus the actual response values and 
detects the values that are not easily predicted by the 
model. The data points were on the 45 degree line that 




^ 1 \ \ \ 

0 5 10 15 20 25 

Actual 

Figure 3 Studentized residuals versus predicted values to 
check the constant error. 

v J 



means the values were not detected. As a result, the 
diagnosis of residuals reveals that there is no statistical 
problem in the model. 

The model presentation 

The model expresses the relationship between responses 
of actual variables and the variables themselves. The 
validated model is presented in actual variables by 
equation (4), 

Y = -221.67 + 211.82;\:i + 174.00;c2 

+ 37.79;\:i;\:2-104.84;v^-86.95;\:^ (4) 

where the actual values of the variables Xi and X2 were 
shown in Table 2 and Y is the alpha. As shown, the 
parameters including linear {xi, X2)> quadratic {xl, xl) 
and interaction {X1X2) affected on the interested res- 
ponse. The effects are presented by the individual co- 
efficients and the related signs (+, -) in the model. The 
coefficients indicate the specific weight of the parame- 
ters in the model while the signs are synergistic (+) and 
antagonistic (-) effects of variables on the response (Y). 

Table 5 The summary of optimized Input variables and 
obtained maximized photodegradatlon% by canonical, 
graphical, numerical methods and validation value of the 
photodegradatlon 



Method 


BiaOs (%mol) 


Ti02 (%mol) 


Alpha 


Canonical (point) 


1.195 


1.025 


15.03 


Grapliical (area) 


Area around 1.25 


Area around 1.25 


20.03 


Numerical (prediction) 


1.24 


1.27 


20.03 


Validated sample 


1.24 


1.27 


21.6 
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Figure 5 The microstructure of the optimized varistor morphology, (a) SEM (b) distribution of ZnO grain size. 



The weights determine the importance of the parame- 
ters roles in the modeling. The model is able to 
optimize input variables and also approximately predict 
the response inside of the actual experimental region 
that confirmed by minimum standard error (Figure 1) 
[33]. Therefore, the model was used to optimize the 
molar ratio of Bi203 and Ti02 to achieve maximized 
alpha. 

The model optimization 

The model is able to optimize the variables by using ca- 
nonical response and graphical plots. The canonical res- 
ponses, local optimums, in terms of the code and actual 



variables were determined by differentiating the model 
(Eq. 4) as presented in equations (5 and 6), 



W/dx,]^^ = 0 



(5) 
(6) 



where the terms were introduced in Table 2. Therefore, 
the optimum canonical amount of Bi203 and Ti02 were 
1.195 and 1.025 respectively. At this optimum, the maxi- 
mized alpha was 15.03. In fact, the optimization is a kind 
of the tradition methods, one-variable-at-a-time, because 
in each case, one of the variables was varied and others 
were constant. In graphical optimization, the model 
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simultaneously considered the effect of Bi203 and Ti02 
on the alpha. Figure 4 presents the three dimension re- 
sponse surface (3D) plot for the synergy between Bi203 
(1-1.5 mol%) and Ti02 (1-1.5 mol%) which is standard 
error limitation area (Figure 1). As shown, the alpha in- 
creased within 1 to 1.25 mol% Bi203 and Ti02. However, 
when the amount of Bi203 and Ti02 was increased in 
excess of the optimum (1.25 mol%), the alpha decreased. 
Therefore, the optimum was determined in a surface 
area around 1.25 mol% for the both additives. The ma- 
ximum alpha was 20.031 at center of the surface which 
indicated by a flag on top of Figure 4. 

The model predicted the maximized alpha in desir- 
able condition of the additives and standard error 
which facilitated by default of software numerical op- 
tion. The desirability is an objective function that uses 
mathematical methods to find the optimum condition 
[25]. The range of the function is from zero (outside 
of the limit area) to one (at the goal). The criteria for 
this case Bi203 (in range), Ti02 (in range), standard 
error (minimized) and alpha (maximized). The sug- 
gested solution was Bi203 (1.24 mol%), Ti02 (1.27 
mol%), standard error (0.35), and alpha 20.03. The de- 
sirability of the solution was 0.981 which is close to 
100% (at the goal). The solution was performed to 
confirm the prediction by validated experiment. The 
validated alpha was determined 21.6 which was quite 
close to the predicted alpha (20.03). Table 5 illustrated 
a summary of optimized molar ratios of Bi203 and 
Ti02 and also the related maximized alpha which 
obtained by canonical, graphical and numerical model 
optimization method. 



The validated varistor 

The validated sample was characterized as final varistor for 
this optimization process. Figure 5 demonstrates the SEM 
morphology of the sintered ceramic microstructure of the 
varistor. Figure 5a indicates the great homogeneity of ZnO 
grain size. The average grain size was 15 (im (Figure 5b). 

Figure 6 illustrates EDX spectra of a limited area of the 
etched varistor surfaces composition. As shown all addi- 
tives particular Bi (1.0 weight%) and Ti (1.17 weight%) 
were detected in the selected area after sintering process 
at 1260°C. 



< 80 




♦ First try 

• Second try 
A Third try 
□ Fourth try 



50 80 110 

E (V/mm) 



140 170 



Figure 8 E-J characteristic curves of the optimized samples for 
first to fourth measurments. 
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Figure 7 shows the XRD pattern of the optimized sam- 
ple which presented three phase of ZnO, spinel and metal 
oxide of additives. The composites were including ZnO 
(00-005-0664), Bi203 (00-002-0988), TiOs (01-072-0020), 
Mn02 (00-003-1041), CoO (00-048-1719), Sb203 (01- 
075-1567), AI2O3 (00-004-0879) and Sb3Ti20io (00-028- 
0103). The number that mentioned in parenthesis are 
XRD reference code. 

The electrical properties of the varistor were basis of 
I-V characteristic measurement that shows breakdown 
voltage was 98 V/mm with non-linear coefficient 21.6. 
The leakage current was 0.013 mA/cm^. The stability 
of the varistor was measured by alpha recovery after re- 
moving the over voltage (Figure 8). As shown the stabi- 
lity was quit significant after fourth over voltage. 

Conclusions 

This work reports modeling and optimization of the molar 
ratio of Bi203 and Ti02 by RSM. The fabrication was 
designed by CCD using two variables and a response. To 
obtain actual responses, the design was performed in 
laboratory by conventional fabrication methods. The ac- 
tual responses were fitted into a quadratic model. The 
model was validated by ANOVA which provided evi- 
dences such as high F-value (153.6), very low p-value 
(<0.0001), Radj (0.985) and Rp^ed (0.947). The results of the 
validation showed the model was significant. The model 
tracked the optimum of the designed additives by using 
3D plots. In the optimum condition, the molars ratio of 
Bi203 and Ti02 were around 1.25 that maximized the 
alpha value at 20. Moreover, the model suggested a solu- 
tion to predict the optimum amount of the additives. In 
this case, the condition of the solution included standard 
error of 0.35, Bi203 of 1.24, Ti02 of 1.27 and alpha of 
20.03. The solution was tested by further experiments. As 
the validation test showed, the obtained value of the alpha 
(21.6) was very close to the predicted value (20.03). There- 
fore, RSM was succeeded in modeling of the additives in 
fabrication of zinc oxide based low voltage varistor to 
achieve maximum alpha. 
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